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Abstract 



We prove the commutativity of the first two nontrivial integrals of motion 

>■ ' for quantum spin chains with elhptic form of the exchange interaction. We 

^ ■ also show their liner independence for the number of spins larger than 4. 

Q> . As a byproduct, we obtained several identities between elliptic Weierstrass 

functions of three and four arguments. 
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1 Introduction 



In this paper we consider the question of integrability of the quantum model 
5^ \ related to the ID Heisenberg chain with non- nearest, variable range exchange 

interaction. It is defined by the Hamiltonian 

TV 

H = ho Y. hij-k)P,k, (1) 

j.fc = 1 
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where A^ > 2 is an arbitrary integer, the transposition operators {Pjk} form 
an arbitrary representation of the permutation group S^, in particular, they 
obey the relations 



jkl ^Ijk 



Pklj, Pjkl = PjkPkl, (2) 



ioT j ^ k ^ I ^ j, ho is a. coupling constant and h{j — k) = p{j — k), where 
p{x) is the elliptic Weierstrass function with real period N and complex 
period u = ik, k E H, k being a free parameter. The model reduces to the 
Heisenberg spin-| chain if 

Pjk = ^{l + ^j-^k), (3) 

where {(Tj} are usual Pauli matrices. 

There is almost no doubt about the integrability of the above Heisenberg 
chain introduced first in [Ij. Indeed, in [1] one of us found the Lax represen- 
tation of the Heisenberg equations of motion for (1) with (2) and the first 
two nontrivial integrals of motion which can be written in compact form 

N 

J= E f{j-k)fik-l)fil-j)P^ki, (4) 



3,k,l = 1 



where 



/(.) = '-^e-^'-^ (5) 

a(x)(T[a] 



and ({x),a{x) are the Weierstrass functions (e.g. [2]) related to p{x) as 

dC{x) d\oga{x) . 

—i— = -p{x), -; = C{x). (6) 

Due to the arbitrariness of the "spectral" parameter a, (4) in fact contains 
only three independent integrals of motion, 

J = -^p'{(^)Jo + p{a)Ji - 2'^2, 

N 

Jo = E -^i^^' 

j, fc, i = 1 



N 



Ji= H ^jkiPjki, 'fjki = C(j -k) + C{k -l) + C(^ - j), (7) 



j,k,l ^ 1 



N 



■h — 2^ FjkiPjki, (8) 



j, fc, i = 1 

3 7^kj^l^3 



Fjki = l{2[aj-k)+ak-i)+ai-j)] 

x[p{j -k) + p{k -l) + p{l-j)] + p'U -k) + p\k -l) + p'{l-j)}. (9) 

The integral of motion Jq, related to the total spin for the case of spin chain 
(3), trivially commutes with any transposition Pj}. and therefore also with 
if, J I and J 2- We are concerned mainly with Ji and J 2- It is easy to prove 
the identities 

2[a3~k)+ak-l)+ai-3)M3-k) + p'U-k)=2[a3-k)+ak-l)+ai-3)]^ 

p{k-i)+p'{k-i) = 2[aj-k)+ak-i)+ai-j)]p{i-j)+p'{i-j), (lo) 

which allow one to rewrite Fjki in one of the following forms: 

F,,i = 2[C(j -k)+ C{k - /) + C(/ - 3)M3 -k) + p'{3 ~k)= (11) 

= 2[C(j -k)+ ({k -l) + C(/ - 3)Ml - 3) + P\l - 3)- (12) 

Note that both ipjki and Fjki are antisymmetric with respect to permutations 
of their indices. 

Unfortunately, the Lax pair formalism cannot produce higher integrals 
of motion due to quantum nature of the problem. The eigenvectors of (1) 
with (2) were explicitly found in [3J up to the solutions of the transcendental 
Bethe ansatz-like equations. In the trigonometric limit of the Weierstrass 
functions (k — >• cx)), one recovers the Haldane-Shastry model [1], and Ji 
might be reduced (for the spin representation (3)) to the product of the 
Yangian generator Y2 = J2f=i,j^k^'^^ 'kU ~ k){aj x cxk) and the total spin 
•S* = I J2jLi <?j iS] • Hence in this limit the symmetry of the model is the 
Yangian Y{sl{2)), and the mutual commutativity of Ji and J2 has been 
proved rather easily [5]. As for general elliptic case, it is highly nontrivial 
problem, and we would like to solve it in this paper. If the commutativity 
does not take place, there would be a whole series of the nontrivial operators 



[Ji, J2], [Ji, [Ji, J2]] etc., commuting with the Hamiltonian (1) (as in the case 
of the components of Y2 which do not commute). Till now, there is no way 
to include the elliptic model (1) into the general quantum inverse scattering 
method |6J. Therefore we shall use the direct method of the evaluation of 
the commutator. 

2 Commutativity of Ji and J2 

Let us write down the commutator of the operators (7), (8) in the form 

N N 

[Jl, J2] = 2^ 2^ V^jklFmnp[Pjkl, Pmnp], (13) 

J,fc,Z — 1 m. n. p ^ 1 

J#fc7^'#J m ^ n ^ p ^ m 

The commutator at the right-hand side of (13) might be nonzero if and 
only if one or two indices {mnp) coincide with (jkl). Consider first the 
coincidence of one index (say, m) with one of (jkl). The direct calculation 
of this contribution to the commutator can be written as 

N 
•^3 = 9 2^ {^jnpFjkl — fjklFjnp)Pjklnp, (14) 

j, k, I, n, p — 1 
all different 

where Pjkinp = PjkPkiPinPnp is symmetric with respect to all cyclic permuta- 
tions of its indices. Hence the coefficient in front of it can be rewritten due 
to this symmetry, and one finds 

9 ^ 

^3 ^ T 2.^ ^f'jklnp^jklnpi \^'^) 

j, k, l^ n^ p ^ 1 
all different 

where 

pjkinp = Fjkl{^jnp — Vlnp) + Fjnpi'^nkl — ^jkl) + Fkin{(pjkp — '^jnp) 

+Fjkp{^pln — '^kln) + FinpiifjM " ^jkp)- (16) 

The function Qjkinp in fact depends on four arguments due to the fact that 
if and F depend only on differences of their indices. Let us introduce the 
notation 

p — j = x,p — k = y,p — l = z,n — p = v. (17) 



Then other differences can be written as 

n—j = v+x,n—k = v+y,n—l = v+z,j — k = y—x,j—l = z—x, k—l = z~y. 

(18) 
With the use of (7), (12), (17-18), we can rewrite Q (16) as 

^jkinp = R{x,y,z,v), 

where 

R{x,y,z,v) = 

[2{C{y-x)+C{z-y)+C{x-z))p{x-z)+p'{x-z)][-C{v+x)+C{x)+C{v+z)-C{z)] 

+[2(C(x) + ({y -x)- ay)My) - p\y)m^) + C(^^) - C(^ -y)- C(^ + y)] 

^2{av)H{x)~av+x))p{v+x)-p\v+x)][av+y)-av+z)-ay-x)+az-x)] 
+[2{-av+z)+av)+az)Mz)+p\z)][ax-z)+az-y)~a^)H{y)] m 

^2{az-y)-av+z)+av+y))p{v+y)+p\v+y)][ay-x)-ay)-av)+av+x)] 

Our goal is now to simphfy this very cumbersome formula. First, let us note 
that R{x, y, z, v) is elliptic, i.e. double periodic function of all its arguments. 
And second, we shall use the following Laurent decomposition of C^{x) and 
p{x) near x=0, the only singularity point of them, 

p{x)^x-^ + ax'^ + 0{x^), C{x)--x~^ --x^ + 0{x^), (20) 

and the differential equations for the Weierstrass p function, 

p\xf = Apixf - g2p{x) - gs, p"{x) = Qp{x)^ - -j, (21) 

where a = ^,g2, gs are some constants. 

Consider now R{x, y, z, v) as the elliptic function of v. It can have simple 
poles at four points: v = 0,v = —x, v = —y, v = —z and no other singulari- 
ties on the torus T = C/{ZN + Zuj). It might be equal to zero if we would 
prove that all these poles are in fact absent (in this case R does not depend 
on v) , and that R{x, y,z,0) = 0. 

Let us calculate the Laurent decomposition of R near the point v = 0. It 
reads 

R{x, y, z, v) ~ v~'A{x, y, z) + B{x, y, z) + ..., (22) 



where 

A{x, y, z) = 2[p(x)(C(t/) - CW - C{y -x) + az - x)) 

+p{z){ax -z)+ az -y)- ax) + ay)) 
+pma^) + cw - 2C(y) + ay -^)+ ay - ^)) - p'm, (23) 

B{x, y, z) = -p'{x){az) - ((y) + ({y - x) - ({z - x)) 

+p'{z){ay) - a^) + a^ -z)+ a^ - y)) 
+p'{y){ax) + cw - 2C(y) - az -y)+ ay - ^)) 

+2{p{x) - p{y)){p{z) - p{y)) - p"{y). (24) 

Consider A{x, y, z) as the elhptic function of the argument x. It might have 
poles ai X = 0,x = y,x = z. Let us calculate the first two terms of its 
Laurent expansion near x = 0: 

A{x, y, z) ~ 2{x-^[x{p{z) - p{y)) + y (p'(t/) - p'{z))] + p{z){-x-' + C(l/) 

+C{z -y)- az)) + p{y){x-' + C(^) - C{y) + C{y - z)) - p'{y)} (25) 

= -ip'iy) + p'{z)) + 2{p{y) - p{z)){az) - ay) + Ciy " z)). 

Now we see that A{x,y,z) has no pole at a; = and A{0,y,z) = due to 
the known identity (change z to —z in the composition formula for () 

lip'iy) + P'{z)) = {p{y) - p{z)){az) - ay) + Ciy - z)). (26). 

An easy calculation based on (23) shows that there are no poles of A at 
X = y and x = z. We conclude that 

A{x,y,z)=0. (27) 

Let us now simphfy the expression (24) for B{x,y,z). It is easy to see that 
aX X = y and x = z there are no poles of this function. The calculation shows 
that there is no pole at x = too and gives 

5(0, y, z) = i(p"(y) - p"iz)) + ip'iz) - p'iy))iay) - Ciz) + az - y)) 

-2p{y){p{z)-p{y))-p"{y). (28) 



By using the identity (26) and differential equations (21) one can write 
-8(0, y, z) in the form 

B{0,y,z) = 2{p{yf - p{zf) - 2p{y){p{z) - p{y)) 

2 p{z) - p{y) 2 

Hence the elhptic function B{x,y,z) has no poles and B{0,y,z) = 0. It 
results in the identity 

B{x,y,z) = 0. (29) 

Let us summarize these steps of calculations. We proved that R{x, y, z, v) 
has no pole at f = and R{x,y,z,0) = 0. But it might have poles at 
V = —X, —y, —z. Calculation of the asymptotics at f -^ ~x gives 

R{x, y, z, v) ~ — -^[2(C(^ -y) + C{x -z) + ({y - x)){p{y - x) - p{x - z)) 

+p'iy-x) + p'iz-x)]. (30) 

But the identity (26) shows that the right-hand side of (30) is just zero. 
Similar calculations result in the absence of poles of R{x, y, z, v) ai v = —y 
and V = —z. Hence this function has no poles in v at all and R{x, y, z, 0) = 0. 
We are coming up to the identity 

R{x,y, z,v) = 0. (31) 

It means that all contributions to the commutator [Ji, J2] quartic in permu- 
tation operators (14) disappear. Let us consider now the case of coinciding 
two pairs of indices in the sets (jkl), {mnp) in (13). The corresponding 
contribution to the commutator consists of two parts, 

N 
Ji = 9 2^ {^jklFjkn — '^jknFjkl)PjlPkn, (32) 

j, k, I, n ^ 1 
all different 

N 

<^5 = 9 2^ Fjkp{ifijp — (pkip)Pjki- (33) 

j, k. I. p — 1 
all different 

The operator in (32) is invariant under changing indices (j ^^ i)',{k ^^ n); 
{j ^-^ k,l 'r^ n). Symmetrization of the coefficient in front of it gives, after 



an easy calculation taking into account the antisymmetry of ^pjki and Fj^i 
under the transposition of two indices, that J4 = for otherwise arbitrary 

fjM, Fj ki- 
lt remains to calculate J5. Since Pjki is symmetric with respect to the 
cyclic permutations of (jkl), (33) can be written in the form 

N 

"^5^3 / ^ -L jklp^jkh 

j, k,l, p ^ 1 
all different 

Tjkip = Fjkp{(pijp — <Pkip) + Fi,jp{ifkip — (pjkp) + Fkip{(pjkp — ^ijp)- (34) 
Let us introduce the notation 

j— p = x,k— p = y,l— p = z. (35) 

Then 

j — k = x — y,k — l = y — z,j— l = x — z, (36) 

and we can rewrite (34) with the use of (7), (11) as 

Tjkip = ^ix,y,z) = (37) 

[2(C(x -y) + C{y) - C{x))p{x - y) + p'{x - y)] 

X [({z ~x) + ({z -y) + C(x) + Ciy) - 2C(^)] 
+ [2(C(2; - x) + C,{x) - C,{z))p{z - x) + p\z - x)\ 

X Uy - ^) + C(y - :^) + C(a:) + C(^) - K{y)\ 
+[2(C(y - z) + C(^) - C(i/))p(y - ^) + p'(i/ - ^)] 

X [C(a: - y) + C(:r - ^) + C(y) + C(^) - 2C(a:)]. 

It is easy to see that $(x, y, z) is antisymmetric with respect to permutations 
of its arguments, 

$(3;, y, z) = -$(1/, X, z) = -$(x, z, y). (38) 

The problem consists now in simplifying ^[x^y^z^ which is elliptic function 
of all its arguments. As a function of x, it has poles at a; = 0, x = y, a; = z. 
Let us calculate the first two terms of its Laurent expansion near x = 0, 

$(x, y, z) ~ 2x-2[p(^) - p[y)] + x-^{2[az - y) - C(^) + C(l/)] 



x[2p{y -z)- p{y) - p{z)] + p'{y) - p'{z) - 2p'{y - z)}. 

The coefficient at x^^ can be drastically simplified by using the identity (26). 
Implying it two times results in 

^x,y,z) ~ 2x-'[piz) - piy)] + 2x-'[p'iy) - p'iz)]. (39) 

The first coefficients in the Laurent expansions near the points x = y and 

X = z are 

<^{x,y,z)r^-2{x-y)-^p'{y), <^{x,y,z) r^ 2{x - z)-^p'{z). (40) 

Let us consider now the trial function 

^{x,y,z) = 2{p{x - y)[p{y) - p{x)] + p{x - z)[p{x) ~ p{z)]}. (41) 

It is easy to see that it has poles at x = 0,x = y,x = z with the same 
residues (39), (40) as ^{x,y,z). Hence 

^{x, y, z) = ^(x, y, z) + tpiy, z), 

where ilj{y,z) does not depend on x. Now, with the use of antisymmetry of 
$ (38), one finds that the only choice for ip is 

i^{y, z) = 2p{y - z)[p{z) - p{y)] 

and finally one can write the remarkable identity 

^{x, y, z) = 2{p{x - y) [p{y) - p{x)\ 

+p{x - z)[p{x) - p{z)] + p{y - z)[p{z) - p{y)]}. (42) 

Now let us prove the relation 

N 

J2 Hx,y,z) = (43) 

for any fixed j ^ k ^ I ^ j. Indeed, coming back to the notation (35-36) 
and using (42), one finds 

N 

E Hx,y,z) = 2{p{j-k)[q{k,l,j)-q{j,k,l)]+p{j-lMj,k,l)-q{l,j,k)] 



+p{l~k)[q{l,j,k) -q{k,l,j)], 
where 

N 

<l{kJJ)= J2 p{k~P) = S{k)-p{k-j)-p{k-l), 
p = 1 

P 7^ j> *:> i 

N 



s{k)= y: pik-p). 



I = 1 

^ k 



But S{k) does not depend on k since p{k — p) is periodic with the period A^. 
Now it is easy to see that (43) holds for all j, k, I and the commutator [Ji, J2] 
vanishes. 

3 Linear independence of Jq, Ji and J2 

Let us prove now that the integrals of motion Jq, Ji and J2 are linearly 
independent for A^ > 4. More specifically, we prove that the operator Jq is 
linearly independent of Ji and J2 for A^ > 3, operators Ji and J2 are linearly 
dependent for A^ = 3,4, and operators Ji, J2 are linearly independent for 

N >4:. 

To study the linear independence, we are looking for the complex numbers 
A, fj,, p such that 

A Jo + /iJi + pJ2 = 0. 

As the coefficients in equations (7) and (8) are symmetrized with respect to 
the cyclic permutations of indices, the last relation is equivalent to 

A + pi-Pjki + pFjki = 

for any mutually different j,k,l = 1, . . . ,N. As ipjki and Fj^i are antisym- 
metric under the exchange of two indices, this is further equivalent to 

A = , fupjki + pFjki = 0. 

In particular, Jq is linearly independent of Ji and J2. 
Let us now consider the case of A^ = 3. Here 

Jl = 3(/9i23(Jl23 ~ J213) 5 J2 = 3-Fl23(Jl23 ^ J213) 

10 



so Ji and J2 are linearly dependent for A^ = 3. In the case of A^ = 4, we 
obtain remembering that A^ is the period of Weierstrass functions in our 
considerations and their other properties 

Ji = 3v9i23(-Pi23 — -P213 + -P124 — -P214 + -P134 — -P314 + -P234 — -P324), 

J2 = 3Fi23(Pl23 — -P213 + -P124 — -P214 + -P134 — -P314 + -P234 — -P324) 

with 

99123 = C(2) - 2C(1) , Fi23 = ^[2y.?23-P'(l)] 

and the linear dependence of Ji and J2 is seen for A^ = 4. 

Let us further on assume A^ > 4 and assume that there exists /j, and p 
satisfying equations fupjki + pFjki = for every possible j, k, I. Let us fix k 
and / and define a functions 

^(z) = pa{z) + pb{z), 

a{z) = C{z -k) + C{k -l) + C(/ - z), 
c{z) = p'{z -k) + p{k -l) + p'{l - z), 



K^) = I ' 



c(z) + 2a(zf 



such that our equations read ip{j) = for j G {1, . . . , A^} \ {k, I}, ip is an 
elliptic function with periods A^ and u. The only possible poles are at the 
points z = k and z = I. They are simple poles of a, let us calculate the 
behavior of h{z) for z = fc + x, x — > 0: 

a{z) = -- C'(/ - k)x + 0{x^) = - + p{l- k)x + 0(x2), 

X X 

2 
c{z) = — ^ + 0(x), 

X'^ 

2p([-fc)^ 

X 

Similar formulas hold for z -^ I due to the antisymmetry of a, c, h with respect 
to the interchange of k and /. Therefore if) has at most simple poles at k and 
/. By Liouville theorem (e.g. [2]), ^/^ can have at most two zeroes (modulo 
periods) if it is not constant. However, there are at least N — 2 > 2 zeroes 



11 



at the points z = j. So ip{z) = 0. Looking for the behavior at z ^ k, we see 
that equation 

H + 2p{l -k)p = 

must be vahd. The function p can take the same value at most twice (modulo 
periods) due to the Liouville theorems. Ask ^ I can be chosen arbitrarily, we 
find two different values among N —1 numbers p{l — k) for l — k = 1, . . . , A^— 1. 
So necessarily p = p = and the hnear independence of Jq, Ji, J2 is proved 
for A^ > 4. Their linear independence of iJ is trivial as different permutations 
enters the definition of H. 

4 Conclusions 

To summarize, we proved that the Hamiltonian (1) and operators Jq, Ji (7) 
and J2 (8) are linearly independent and generate the commutative ring. As a 
byproduct, we obtained the remarkable identities between elliptic functions 
(31), (42). The proof was based on direct evaluation of [Ji, J2] due to the 
lack of any other methods. The model (1) with elliptic form of h{j — k) is 
still not immersed in the scheme of the quantum inverse scattering method. 
This is highly desirable task which we postpone for further study. The pres- 
ence of the operators of higher orders in permutations commuting with the 
Hamiltonian was also mentioned [7] but till now there is no way to prove 
their mutual commutativity. 
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